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Abstract 

Background: DNA methyltransferase 3A (DNMT3A) mutations were considered to be independently associated with 
unfavorable prognosis in adults with de novo acute myeloid leukemia (AIVIL), however, there are still debates on this topic. 
Here, we aim to further investigate the association between DNMTBA mutations and prognosis of patients with AIVIL. 

Methods: Eligible studies were identified from several data bases including PubMed, Embase, Web of Science, ClinicalTrials 
and the Cochrane Library (up to June 2013). The primary endpoint was overall survival (OS), while relapse-free survival (RFS) 
and event-free survival (EPS) were chosen as secondary endpoints. If possible, we would pool estimate effects (hazard ratio 
[HR] with 95% confidence interval[CI]) of outcomes in random and fixed effects models respectively. 

Results: That twelve cohort studies with 6377 patients exploring the potential significance of DNMT3A mutations on 
prognosis were included. Patients with DNMT3A mutations had slightly shorter OS (HR= 1.60; 95% CI, 1.31-1.95; P<0.001), 
as compared to wild-type carriers. Among the patients younger than 60 years of age, DNI\/1T3A mutations predicted a worse 
OS (HR = 1.84;95% CI, 1.36-2.50; P<0.001). In addition, mutant DNIV1T3A predicted inferior OS (HR = 2.30; 95% CI, 1.78-2.97; 
P = 0.862) in patients with unfavorable genotype abnormalities. Similar results were also found in some other subgroups. 
However, no significant prognostic value was found on OS (HR=1.40; 95% CI, 0.98-1.99; P = 0.798) in the favorable 
genotype subgroup. Similar results were found on RFS and EPS under different conditions. 

Conclusions: DNMT3A mutations have slightly but significantly poor prognostic impact on OS, RPS and EPS of adults with 
de novo AML in total population and some specific subgroups. 
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Introduction 

Acute myeloid leukemia (AML) are a group of heterogeneous 
diseases with respect to biological and clinical outcomes, which 
have been considered to be related to cytogenetic and gene lesions 
in hematopoietic stem or progenitor cells [1,2,3]. It is extremely 
important for patients with de novo AML to assess risk status 
based on cytogenetic and molecular abnormalities. Cytogenetic 
abnormalities have been utilized widely to evaluate the risk status 
of patients [3,4,5,6]. However, only several molecular lesions were 
applied to risk stratification in clinical practice, such as FMS-like 
tyrosine kinase 3 (FLT3) [7,8], Nucleophosminl (NPMl) [9,10], 
CCAAT/enhancer-binding protein alpha (CEBPA) [11,12]. 
Whereas DNMT3A [13,14,15], Isocitrate dehydrogenase 1 
(IDHl), IDH2 [16,17] and TET2 [18,19] have not been fuUy 
assessed. DNMT3A plays a role in de novo methylation of specific 



CpG islands in DNA, and DNMT3A mutations are implicated in 
the pathogenesis and prognosis of patients with AML [20]. 
Recurrent abnormities of the gene include missense, frameshift, 
nonsense, splice-site mutations and part deletion [13] as well. The 
most common type of missense mutations is predicted to affect the 
amino acid at R882, which accounts for more than half the 
population with DNMT3A mutations [1 3,2 1,22]. These mutations 
are highly enriched in the group of adults with an intermediate- 
risk cytogenetic profile or with cytogenetic normal AML (CN- 
AML) [13,22,23,24], but they are absent or scarce in patients with 
a favorable-risk cytogenetic profile [21,23,24]. DNMT3A muta- 
tions having negative impact on OS of patients with AML were 
reported in several studies but not in the others [20,22,25,26], 
Similar discrepancies were observed in terms of RFS and EFS. 
Thus it is necessary to perform a systematic review and meta- 
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Figure 1. PRISMA flow diagram for study review and inclusion. 
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analysis to further clarify the prognostic values of mutant 
DNMT3A in patients with de novo AML. 

Materials and Methods 

1 Literature search and search strategy 

Literature search was conducted in several electronic databases 
including PubMed, Embase, Web of Science, ClinicalTrials and 
the Cochrane Library (up to June 26, 2013). We used the following 
search terms (medical subject headings [MeSH] or key words) 
limited in titles and abstracts: "acute myeloid leukemia", "acute 
myeloblastic leukemia", "acute myelocytic leukemia", "acute 
granulocytic leukemia", "Leukemia, Myeloid, Acute", "AML", 
and "DNA methyltransferase 3A", "DNA (cytosine-5)-methyl- 
transferase 3A", "DNMT3A". No other search terms were added 
in order to obtain more related studies. 

2 Selection criteria 

If meeting the following selection criteria, studies should be 
included: (1) studies published in English; (2) Cohort studies or 
meeting abstracts; (3) studies focusing on prognostic impact of 
DNMT3A mutations in adults with de novo AML; (4) data from 
time-to-event analysis in every study; (5) the latest version of 
multiple publications or overlap of some context from the same 
cohort should be selected, the older version should be used to 
clarify methodology or characteristics of patients. 



Two investigators (R.-x.T. and T.-s.Z.) were asked to identify 
the eligible studies independently. Endnote x6 were used for the 
management of the articles and for removing most of the 
duplicates. Unrelated articles and remaining duplicates were 
excluded by reading the abstracts carefully. The remaining full- 
text articles were retrieved and reviewed carefully to identify the 
eligible studies. Any divergency was resolved by asking a third 
author's opinion (H.-r.F.). The selection process was documented 
in a flow chart recommended in the PRISMA statement. 

3 Risk of bias from included studies 

The quality of eligible studies was assessed independendy by two 
reviewers (R.-x.T. and T.-s.Z.) using the Newcastle-Ottawa-Scale 
(NOS) [27] for cohort studies. The NOS assigns a maximum of 4 
points for selection, 2 points for comparability, and 3 points for 
exposure or outcome. According to the NOS, We classified the 
qualities of these studies into three groups: high (7-9 points), 
intermediate (4-6 points) and low (1-3 points) qualities [27]. 

4 Data extraction 

Two independent researchers (R.-x.T. and T.-s.Z.) reviewed all 
of the articles that met the inclusion criteria. We selected OS as the 
primary endpoints, RFS and EES as the secondary endpoints 
[6,28]. We extracted characteristics of the studies, the correspond- 
ing HRs with 95% CIs from COX multivariable models. 
Unpublished data in some studies were obtained by contacting 
the authors. If the HRs with 95% CIs were not presented in some 
articles, we would extract the data from corresponding Kaplan- 
Meier curves using the methods described by Parmar and Tierney 
[29,30]. 

5 Data synthesis and statistical analysis 

We pooled the HRs and 95% CIs for OS, RFS and EES 
respectively in total population and subgroups, using random- 
effect and fixed-eflFect models simultaneously [31], but only tiie 
estimates from random effects model was selected as the basis of 
our conclusion, because this approach could provide a more 
conservative assessment of the average effect size [32]. In the forest 
plots, pooled HR values>l represented a positive association 
between mutant DNMT3A and prognosis of patients; while 
pooled HR values< 1 represented a negative association. 

6 Assessment of heterogeneity 

We assessed statistical heterogeneity as to the pooled HR for OS 
in total population by visual inspection of the forests plots and by a 
formal statistical test using Chi-square test with a significance level 
at P<0. 1. I^ statistics was also used to quantitatively assess the 
possible heterogeneity (I^>30% represents moderate heterogene- 
ity, I^>75"/() represents considerable heterogeneity) [33]. 

7 Subgroup analysis 

We explored the possible causes of heterogeneity by subgroups 
analysis. We examined the following as important subgroups: (1) 
cytogenetic intermediate-risk subgroup; (2) cytogenetically normal 
subgroup; (3) patients younger than 60 years of age; (4) patients 
older than 60 years of age. (5) favorable risk genotype subgroup, 
namely low risk molecular group, which referred to the patients 
with mutant NPMl and wild-type ELT3-ITD (NPM1+/ELT3- 
ITD— ); (6) unfavorable risk genotype subgroup, also known as 
high risk molecular group, which included the patients (NPM 1 — / 
ELT3-ITD+, NPM1-ELT3-ITD- or NPM1+/ELT3-ITD+); (7) 
patients with DNMT3A mutations at R882; and (8) patients with 
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study ID 



Mutant DNMT3A as a Prognostic Factor in AML 
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Figure 2. Forest plot of the HR for overall survival of all AML patients. DNMT3A mutations versus wild-type DNIVIT3A. l-V Overall: the pooled 
HR with 95% CI using a fixed effects model; D+L Overall: the pooled HR with 95% CI using a random effects model. 
doi:1 0.1 371/journal.pone.0093353.g002 



DNMT.SA mutations at Non-R882. Of note, NPM1/FLT3-ITD 
risk genotype groups were all belong to CN-AML. 

8 Assessment of reporting biases and sensitivity analysis 

Sensitivity analysis was performed to detect the robustness of 
pooled HR for OS of all patients. Firstly, to test for publication 
bias, We selected OS of total population to form a funnel plot 
[34,35] and contour enhanced funnel plot [36], and undertook the 
Egger test as well. Different types of funnel plots were presented in 
a similar form of scatter plots of study effect sizes (x axis) against 
estimated standard errors (y axis). Contour enhanced funnel plots 
display areas of statistical significance on a funnel plot. Contour 
lines represent diffc-rc-nt Icvch of statistical significance (e. g. 5%, 
10%). If .studies seemed to be missing in areas of non-significance 
(P>10%), the asymmetry of the funnel plot might be due to 
publication bias, although other explanations should still be 
considered. Meanwhile, we attempted to find other reporting 
biases such as outcome reporting bias, language bias, citation bias 
and location bias and so on. 

In addition, we tried to figure out "missing" studies and assessed 
the robustness of pooled HR for OS of the overall population by 
combining contour enhanced funnel plots with the trim and fill 
adjustment method(random and fixed effects linear estimator) 
[37,38,39], All statistical analyses were performed by the software 
STATA (version 12. 0). 



Results 

1 Search results 

As shown in Figure 1, 462 records were obtained by searching 
the databases. After exclusion of 154 duplicates, 308 records were 
further screened by reading their titles and abstracts, and then 241 
records were excluded, in which DNMT3A mutations were not 
regarded as the primary research problems. The full texts of the 
remaining 67 studies were read carefully, and 55 of them including 
32 non-cohort studies, 20 unrelated meeting abstracts and 3 
studies having no sufficient data were further excluded. Eventu- 
ally, 12 studies met the predefined selection criteria, which 
included 2 meeting abstract [40] and 10 original articles 
[1 3,20,2 1 ,22,23,24,25,41 ,42,43] . 

2 Risk of bias in the included studies 

We assessed risk of bias (see methodological quality domains), 
according to nine items of NOS. The qualities of 9 (75%) studies 
were considered high, and the remaining 3 (25%) studies were 
thought moderate. More details were shown in Table 1 and 
Table 2. 

3 Characteristics of the included studies 

The main characteristics of the patients with or without 
DNMT3A mutations in 12 included studies were presented in 
Table 1 and Table 2 respectively, additional characteristics of 
the patients were shown in Tables SI and S2. Of these studies, 5 
originated from America, 4 from Europe and 3 from Asia. 
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Figure 3. Forest plot of the HRs for relapse-free or event-free survival of AML patients. DNIVITSA mutations versus wild-type DNMT3A. l-V 
Subtotal represented the pooled HRs with 95% CIs using a fixed effects model; D+L Subtotal represented the pooled HR with 95% CI using a random 
effects model. 
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4 Data extracting 

We obtained directly most of data from the included studies 
except for one set of data [25], which were extracted indirectly 
from the corresponding Kaplan-Meier curves related to mutant 
DNMT3A. Additionally, we managed to obtain some unpublished 
data by contacting the authors. Hou et al. olfered the HRs with 
95% CIs for OS and RFS in unfavorable and favorable risk 
genotype subgroups with AML. Fabiana et al. offered the HRs 
with 95% CIs for OS, EFS and RFS in subgroups as mentioned 
above. 

5 Association of DNiVlT3A mutations with pretreatment 
clinical and molecular characteristics 

More patients with DNMT3A mutations were found with 
mutations in NPMl, FLT3-ITD and IDHl/2 
[13,20,21,22,23,24,42], whereas CEBPA mutations hardly coex- 
isted in the patients with DNMT3A mutations [20,21,22,42,43]. 
DNMT3A mutations were mainly enriched in the patients with 
CN-AML or with intermediate-risk AML. Patients with favorable- 
risk cytogenetics had scarcely DNMT3A mutations. Comparing 
with wild-type DNMT3A patients, DNMT3A mutated patients 
were older, with higher median WBC counts and percentage of 
bone marrow blasts at diagnosis. Whereas, platelet counts of 
patients were not considered to be associated with DNMT3A 
mutations in most of the included studies [21,22,23,24]. 



6 Association of DNMTSA mutations with clinical 
outcomes in total population 

As shown in Figure 2, data were extracted from 12 studies, 
with a total of 6377 patients, including 1161 patients with 
DNMT3A mutation. In the total population, the patients with 
DNMT3A mutations had a shorter OS as compared to those with 
wild-type DNMT3A (HR= 1.60, 95% CI, 1.31-1.95). There was 
considerable heterogeneity as shown by testing (I^ — 74.4%, P< 
0. 001). 

As shown in Figure 3, data were extractable from 5 out of 1 2 
studies, with a total of 3291 patients, including 715 DNMT3A 
mutated patients. In this population, patients with DNMT3A 
mutations had a shorter RFS as compared to those with wild-type 
DNMT3A (HR= 1.79, 95% CI, 1.25-2.56). There was consider- 
able heterogeneity by f testing (I^ = 75.9%, P = 0.002). 

As shown in Figure 3, data were extractable from 4 out of 1 2 
studies, with a total of 1736 patients, including 210 DNMT3A 
mutated patients. In this population, patients with DNMT3A 
mutations had a shorter EFS as compared to them with wild-type 
DNMT3A (HR= 1.88, 95% CI, 1.40-2.52). There was moderate 
heterogeneity by I^ testing (1^ = 44.3%, P = 0.146). 

7 Association of DNMTSA mutations with clinical 
outcomes in different subgroups 

We pooled the HRs in frxed- and random-effect models 
simultaneously, which were shown in Table 3 and Table 4 
respectively, corresponding forest plots were presented in Figures 
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SI to S8. In the favorable risk genotype subgroup, mutant 
DNMT3A had no prognostic impact on OS or RFS other than 
EFS. And in the non-R882-DNMT3A mutations subgroup, 
mutant DNMT3A had no prognostic value for OS. However, 
DNMT3A mutations were independently prognostic factors for 
OS, RFS and EFS of patients in the remaining 6 subgroups. There 
was considerable statistical heterogeneity detected by test in 
most of these subgroups. 

8 Publication bias and sensitive analysis 

We utilized the data about OS of total population to search 
publication bias and do sensitive analysis. 

As shown in Figure 4, there was a strong suggestion of 
asymmetry in the funnel plot (Egger's test P<0. 001). 

Figure 5 demonstrated the distribution of 1 2 real studies and 4 
filled studies. Of the 12 real studies, 8 studies lay above the right 
5% significant contour line, and 4 studies located below the 5% 
significant contour line. Only one study distributed on the left- 
hand side of the plot, which also implied asymmetry of the plot. 
Secondly, the trim and fill method imputed a total of 4 "missing" 
studies in a random effects model, which were presented in the 
form of gray triangles located in the region of statistical non- 
significance (P>10%), and the vertical black line showed the 
pooled log HR as to OS in 12 studies, while the vertical gray line 
showed the pooled log HR including the filled studies. After 
adding the "missing" studies to the contour-enhanced funnel plot, 
the pooled fixed-effects HR was 1.282 (95% CI, 1.174-1.400), and 
the pooled random-effects HR was 1.394 (95% CI, 1.161-1.674), 
which implied DNMT3A mutations had a slightiy negative impact 
on OS in total population. 

Discussion 

We have performed a systematic review and meta-analysis in an 
effort to figure out the association between DNMT3A mutations 
and poor prognosis of patients with de novo AML. 

So far, the 12 published studies have addressed the prognostic 
value of mutant DNMT3A in AJvIL patients. According to the 
rules of NOS, the majority of included studies belonged to high 
quality, only 3 of the studies had moderate quality. A total of 6377 
patients of the 12 included studies, where there were 1161 patients 
with mutant DNMT3A. 

The primary end point was OS. DNMT3A mutations had 
slightly prognostic impact on OS in total population, which was 
consistent with the pooled outcomes in the meta-analysis by 
Shivarov et al [44]. And the conclusion was suitable for the 
majority of subgroups except for the favorable risk genotype and 
non-R882 DNMT3A mutation subgroups. The second endpoints 
were RFS and EFS. First of all, mutant DNMT3A had poor 
prognostic impact on RFS in nearly all of the mentioned groups 
other than favorable risk genotype subgroup. As for EFS, 
DNMT3A mutations were associated with poor prognosis of 
patients in all of the mentioned groups above. However, most of 
the included studies showed neither type of DNMT3A mutation 
had an impact on the probability of achieving CR in multivariate 
analyses [21,22,23,24,41,42]. 

There were significant correlation of DNMT3A mutations with 
patient characteristics, including age of patients and WBC in 
peripheral blood at diagnosis, the mutations of NPM 1 , FLT3 or 
other genes (IDHl, IDH2, TET2, WTl and CEBPA), regardless 
of the types of mutation. In addition, only one study mentioned the 
association between DNMT3A mutations and expression of CD 
molecules on leukemic initiating cells, such as CD 14, CD 13 and 
CD34 [24]. 
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Funnel plot with pseudo 95% confidence limits 
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Figure 4. Funnel plots illustrated significant asymmetry on HR for overall survival of all patients. Studies were distributed 

asymmetrically and suggested biases exist. 

doi:10.1371/journal.pone.0093353.g004 



We observed considerable or moderate between-study hetero- 
geneity in total population and subgroups, unlike that in the meta- 
analysis by Shivarov et al [44]. 

The heterogeneity mainly came from one study lately published 
[26]. It was worth noting that we did not conduct a meta- 
regression analysis because of fewer studies included. So we 
explored subjectively the clinical heterogeneity as follows: firstly, 
with respect to FAB classification, patients' years of age, 
cytogenetic and molecular abnormalities, the constituent ratio of 
patients were various in each study, for instance, 7 studies only 



included patients younger than 60 years of age 
[20,21,23,40,41,42], 2 .studies only included cytogenetic normal 
patients [22,42], 1 study included patients belonging to cytogenetic 
intermediate risk group [22]; secondly, most of the included 
studies did not presented the ECOG performance status of 
patients except for two of them [23,43], so we could not further 
assess the status of patients when they were included, which might 
be one of the sources of heterogeneity; thirdly, the treatment 
programs such as induced chemotherapy and hematopoietic stem 
cell transplantation were diverse in individual studies, even 
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Figure 5. Confunnel with filled studies from metatrim: mutant DNIVIT3A versus wild-type DNMT3A in a random effects model. The pooled HR 
on overall survival from 12 published studies is robust and the heterogeneity mainly results from unpublished studies. 
doi:1 0.1 371/journal.pone.0093353.g005 
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intensities of the same induced chemotlierapy drug were various in 
different studies, so the selection of different treatment programs 
might lead to the heterogeneity of clinical outcomes; eventually, 
there were much between-study heterogeneity in other aspects 
such as the time of follow up, percentages of bone marrow blasts, 
races of patients, types of DNMT3A mutations, the methods of 
gene sequencing, and the number of adjusted covariates among 
the 12 studies and so on. 

We found the asymmetry in a usual funnel plot. Subsequently, 
we adjusted pooled HR for OS in total population by combining 
the trim and fill method with contour enhanced funnel plot, and 
found that the asymmetry was probably caused by unpublished 
"missing" studies, so publication bias was a major threat to the 
validity of our estimate outcomes, although other explanations 
should still be considered. Meanwhile, we verified the robustness 
of the pooled random-effect and fixed-effect outcomes of 16 
studies consistent with the pooled outcomes of 1 2 original studies, 
which was extremely important for our eventual conclusion. 

We did not find any methodological issues in the preparation of 
the review which could put it at risk for bias. 

There are the following aspects of main limitation in our review. 
Firstly, the "missing" studies have been relatively absent to date. 
Secondly, selective reporting were seen in some studies, for 
instance, the pooled HRs for RFS and EFS were shown in 
relatively fewer studies as compared to them for OS, and some 
subgroups such as older patients were not analyzed in most studies, 
which resulted in losses of valuable information. Thirdly, two 
studies were published as meeting abstracts, so we could not assess 
its potential risk of bias in detail. Fourthly, the included studies 
were all published in English, so language bias existed. Fifthly, we 
could not access the association between mutant DNMT3A and 
prognosis of Secondary AML because of the lack of efficient data. 
Finally, if these trials could offer individual patient data, our 
analysis would have been much more perfect. 

In conclusion, using some advanced statistical methods and high 
quality evidences, we confirmed objectively the robustness of 
association between mutant DNMT3A and prognosis of adults with 
de novo AML, although there are considerable heterogeneity in 
many conditions. Combining other important genetic biomarkers, 
DNMT3A mutations would contribute to a more precise clirrical 
risk stratification and decision of treatment. More cohort studies 
concerning DNMT3A mutations arc lux-ded in an effort to further 
verify or modify the pooled estimates to a certain extent, especially 
for the patients with CN-AML or with cytogenetic intermediate risk 
abnormalities. It is meaningful to clarify the mechanism that mutant 
DNMT3A function as a epigenetic regulator in adult AML. And 
more studies will be needed to determine whether dose-intensified 
induction chemotherapy improves the survival outcomes in patients 
with DNMT3A mutations. 
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